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Abstract 



The crosscap number of a knot is an invariant describing the non-orientable surface of 
smallest genus that the knot bounds. Unlike knot genus (its orientable counterpart), cross- 

Ocap numbers are difficult to compute and no general algorithm is known. We present three 
methods for computing crosscap number that offer varying trade-offs between precision and 
speed: (i) an algorithm based on Hilbert basis enumeration and (ii) an algorithm based on 
exact integer programming, both of which either compute the solution precisely or reduce it 
to two possible values, and (iii) a fast but limited precision integer programming algorithm 
that bounds the solution from above. 

The first two algorithms advance the theoretical state of the art, but remain intractable 

for practical use. The third algorithm is fast and effective, which we show in a practical 

^> setting by making significant improvements to the current knowledge of crosscap numbers in 

f"\] | knot tables. Our integer programming framework is general, with the potential for further 

00 ■ applications in computational geometry and topology. 

CO 



1 Introduction 



Knot invariants lie at the heart of computational knot theory. Nevertheless, computing invari- 
ants can be challenging: algorithms often require complex implementations and exponential 
running times, and for some invariants no general algorithm is known. 

In this paper we focus on invariants of knots in K 3 that relate to 2-dimensional surfaces: 
knot genus and crosscap number. In essence, these measure the simplest embedded orientable 
and non-orientable surfaces respectively that have a single boundary curve following the knot. 

Knot genus is well-studied: algorithms are known [15] , and precise values have been computed 
for all prime knots with < 12 crossings 10J . In contrast, although crosscap numbers can be 
computed for special classes of knots jTTJ [TBI &7\ , no algorithm is known for computing them 
in general. Of all 2977 non-trivial prime knots with < 12 crossings, only 289 have crosscap 
numbers that are known precisely [10] . 

The crosscap number displays unusual behaviours that knot genus does not, and embodies 
different information [11, 26 . It is therefore desirable to compute crosscap numbers in a general 
setting. Here we develop three algorithms with varying precision-to-speed trade-offs, which yield 
both theoretical and practical advances. 

Our algorithms do not guarantee to compute the crosscap number precisely for every input 
(this remains an open problem), but they do come close. The first two algorithms are significant 
theoretical advances: they either compute the crosscap number precisely or reduce it to one of 



two possible values. The third algorithm is a significant practical achievement: although it only 
outputs an upper bound on the crosscap number, its strong practical performance combined 
with known lower bounds allows us to make significant improvements to crosscap numbers in 
existing tables of knots. 

The first algorithm, described in Section @J uses Haken's normal surface theory [T5j to reduce 
the computation of crosscap number to a Hilbert basis enumeration over a high-dimensional 
polyhedral cone. Although this general approach follows a common template in computational 
topology, there are complications that cause the usual theoretical techniques to fail. To address 
this, we introduce a special class of triangulations called suitable triangulations, described in 
Section [3] with which we are able to solve these theoretical problems. 

The second and third algorithms, described in Sections [5] and [6j draw on techniques from 
discrete optimisation theory. Here we develop an integer programming framework that allows us 
to approach problems in computational topology using off-the-shelf optimisation software. The 
difference between the second and third algorithms is that the second requires expensive exact 
integer arithmetic; the third makes concessions that allow us to use fast off-the-shelf solvers 
based on floating point computation, but with the side-effect that it only outputs an upper 
bound. 

The first two algorithms remain intractable for all but the simplest knots, though each 
has different strengths through which it may become practical with the growth of supporting 
software in algebra and optimisation. In contrast the third algorithm is extremely fast, and in 
Section[7]we run it over the full 12-crossing knot tables. The results are extremely pleasing: of 
the 2688 knots with unknown crossing number, for 747 we can improve the best-known bounds, 
and for a further 27 we can combine our output with known lower bounds to compute the 
crosscap number precisely. 

The integer programming framework that we introduce here is general, and has significant 
potential for use elsewhere in computational geometry and topology. We discuss these matters 
further in Section [7] 

2 Preliminaries 

Here we give a short summary of concepts from knot theory and normal surface theory that 
appear within this paper. This outline is necessarily brief; for details see the excellent overview 
in [H]. 

For our purposes, a knot is a piecewise linear closed curve embedded in M 3 . We also treat 
knots as being embedded in the 3-sphere S 3 , where S 3 is the one-point compactification M 3 U 
{oo}. A knot K is typically given as a knot diagram, which is a general position projection of 



the knot onto the plane as illustrated in Figure 1(a) See [TS] for more precise definitions of 
these concepts. 





(a) (b) 

Figure 1: A knot diagram and a spanning surface for the trefoil knot 



We are interested in properties of knots related to two-dimensional surfaces. In this paper all 
surfaces are piecewise linear, and might be disconnected unless otherwise stated. A surface with 
no boundary curves (such as a sphere or a torus) is called closed, and a surface with boundary 
curves (such as a disc or a Mobius band) is called bounded. 

Every closed orientable surface is topologically a 2-sphere with g > orientable "handles" 
attached (i.e., a g-holed torus); the orientable genus of such a surface is g. Every closed non- 
orientable surface is a 2-sphere with g > 1 "crosscaps" attached; the non- orientable genus of 
such a surface is g. We simply use the word genus when orientability is clear from context. If an 
(orientable or non-orientable) surface S has boundary, then the (orientable or non-orientable) 
genus of S is the genus of the closed surface obtained by filling each boundary curve of S with 
a disc. 

Let V be a polygonal decomposition of a surface S. The Euler characteristic of V is defined 
as V — E — F, where V, E and F represent the number of vertices, edges and faces of V 
respectively. Euler characteristic is a topological invariant of a surface, and is denoted by x(S). 
The Euler characteristic of an orientable surface of genus k with b boundary curves is 2 — 2fc — b, 
and the Euler characteristic of a non-orientable surface of genus k with b boundary curves is 
2-k-b. 

The knot invariants that we study here relate to surfaces embedded in 3-dimensional space, 
which leads us to study 3-manifolds. A 3-manifold is a higher-dimensional analogue of a surface: 
each interior point of a 3-manifold M has a local neighbourhood that is topologically similar 
to R 3 , and each point on the boundary of M has a local neighbourhood that is topologically 
similar to the closed half-space K 3 > - Again, see [TB] for precise details. A surface ScMis 
embedded in M if it has no self-intersections, and properly embedded if in addition the boundary 
of S lies within the boundary of M, and the interior of S lies within the interior of M. 

Consider a knot K embedded in S 3 . A spanning surface for K is a connected embedded 



surface in S 3 whose boundary is precisely K. Figure 1(b) illustrates a spanning surface for the 
trefoil knot (this surface is a Mobius band, with non-orientable genus 1). 

We can now define the following two invariants of a knot K . The genus of K, denoted 
g(K), is the smallest k for which there exists an orientable genus k spanning surface for K. 
Likewise, the crosscap number of K, denoted C(K), is the smallest k for which there exists a 
non-orientable genus k spanning surface for K . As a special case, the crosscap number of the 
trivial knot (also called the unknot) is defined to be 0. The unknot is the only knot with genus 0, 
and the only knot with crosscap number 0. A key relation between genus and crosscap number 
is the following [11) : 

Theorem 2.1 (Clark's inequality). For any knot K, it is true that C(K) < 2g(K) + 1. 

Let K C S 3 be a knot, and let R be a small regular neighbourhood of K in S 3 . The 
complement of K, denoted K, is the closure of S 3 \R. This is a 3-manifold with boundary, 



obtained by "eating away" the knot from S 3 , as illustrated in Figure 2(a) The boundary 
surface of K is a torus, and any curve on this torus that bounds a disc in R is called a meridian. 
We can reformulate spanning surfaces in terms of knot complements. Consider a knot K C S 3 
and some meridian m on the boundary of K. A spanning surface in K is a connected, properly 
embedded surface S C K with precisely one boundary curve c, where c and m have algebraic 



intersection number ±1 on the torus boundary of A"|^| See Figure 2(b) for an illustration. This is 



essentially the same as our previous definition, since any such surface can be extended through 



1 Essentially we require that the boundary of S cuts m precisely once. Using algebraic intersection number 
allows us to account for any extra trivial "wiggles" back and forth across m. 
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(a) (b) 

Figure 2: The complement of the trefoil knot and a spanning surface within it 




R to give a connected embedded surface bounded by K and vice versa, assuming that the 
neighbourhood R is sufficiently small. 

All of our algorithms work with triangulations of the complement K. In this paper, a 
triangulation of a 3-manifold is a collection of n tetrahedra, some of whose An faces are affinely 
identified in pairs. This broad definition allows for smaller triangulations than a traditional 
simplicial complex; in particular, most triangulations in this paper are one-vertex triangulations, 
where all An tetrahedron vertices are (as a result of the face gluings) identified to a single point 
in K. 

More generally, if T is a 3-manifold triangulation, all tetrahedron vertices that are identified 
to a single point in T are collectively referred to as a single vertex of T; similarly for edges and 
faces. Any vertex, edge or face that lies in the boundary of the underlying 3-manifold is called 
a boundary vertex, edge or face; all others are referred to as internal. Note that the boundary 
faces of T are precisely those tetrahedron faces that are not paired with some partner face. 

In our algorithms we describe spanning surfaces using normal surface theory. A normal 
surface in T is a properly embedded surface that meets each tetrahedron A of T in a disjoint 
collection of triangles and quadrilaterals, each running between distinct edges of A, as illustrated 
in Figure [3] There are four triangle types and three quadrilateral types according to which edges 
they meet. Within each tetrahedron there may be several triangles or quadrilaterals of any given 
type; collectively these are referred to as normal discs. 




Figure 3: Normal triangles and quadrilaterals within a tetrahedron 
The vector representation of a normal surface S is the 7n-dimcnsional integer vector 

V(S) = (fl,l,il,2, tl, 3) *l,4i 91,1191,27 91,3; fe,l, £2,2, £2, 3, ^2,4, 92,1,92,2,92,3; ■ ■ • ,9n,3 ) G Z , 

where tij is the number of triangles in the ith tetrahedron of the jth type, and 9^ is the number 
of quadrilaterals in the ith tetrahedron of the kth type (1 < i < n, 1 < j < 4, 1 < k < 3). 

Theorem 2.2 (Haken, 1961). A vector v = (ti,i, . . . , q n ,3) S K 7 ™ is the vector representation of 
a normal surface if and only if: (i) all elements of v are non-negative integers; (ii) v satisfies 



a certain set of homogeneous linear equations derived from the triangulation, and (in) for each 
i, at most one of the three quadrilateral coordinates <7i,i, gi,2, <7i,3 is non-zero. 

The homogeneous linear equations in (ii) are called the matching equations, and condition (iii) 
is called the quadrilateral constraints. The set of all points in R 7 " with non-negative coordinates 
that satisfy the matching equations is called the normal surface solution cone. 

Given a vector that satisfies all of the constraints of Thcorcm l2.21 the corresponding normal 
surface can be reconstructed uniquely (up to normal isotopy). If X and Y are both normal 
surfaces in some triangulation, the normal sum X + Y is the normal surface with vector rep- 
resentation v(X) + v(Y). It is possible that v(X) + v(Y) does not satisfy the quadrilateral 
constraints, in which case (for our purposes) the sum X + Y is not defined. 

A fundamental normal surface S is one that cannot be written as S = X + Y for non-empty 
normal surfaces X, Y. There are finitely many fundamental normal surfaces in a triangulation, 
corresponding precisely to the vectors in the Hilbert basis of the normal surface solution cone. 

3 Suitable triangulations 

To overcome theoretical difficulties that arise with non-orientable spanning surfaces, we intro- 
duce a special class of triangulations for our algorithms to use. 

Definition. A suitable triangulation T of a knot complement is one for which: 

(i) T has precisely one vertex (and therefore the torus boundary of T contains this one vertex, 
three edges and two faces) p 

(ii) one of the boundary edges of T is a meridian. 

Given a triangulation 7" of a knot complement, it is easy to test whether T is suitable. 
Condition (JTJ) can be verified by grouping vertices of tctrahcdra into equivalence classes under 
identification. For condition (|u|l. we can verify that a boundary edge e is a meridian by attaching 
a solid torus in an appropriate fashion and testing whether the resulting closed manifold is a 
3-sphere|f| 

Our first two algorithms for computing crosscap number require a further condition: 

Definition. An efficient suitable triangulation T of a knot complement is a suitable triangula- 
tion that contains no embedded normal 2-spheres. 

This extra 2-sphere condition is related to, though weaker than, the O-efficiency criterion 
of Jaco and Rubinstein [21] . Testing for efficient suitability is possible but slow: to verify the 
absence of embedded normal 2-spheres one must typically enumerate all extreme rays of the 
high-dimensional normal surface solution cone |21j . 

We give two algorithms for producing a suitable triangulation of a knot complement. The 
first is general but slow; the second is heuristic in nature but works extremely well in practice. 

Algorithm 3.1. Given a knot diagram describing the knot K C S 3 , the following procedure will 
output an efficient suitable triangulation of K . 



2 This follows by a standard Euler characteristic argument. 

3 The solid torus must be attached so that its meridional disc is bounded by e; this can be done using a (1, 1, 0) 
layered solid torus as described in |22| . Although the subsequent 3-sphere test requires worst-case exponential 
time, it is found to run surprisingly fast in practice when the right algorithms and simplification heuristics are 

used nn 



1. Test whether K is the unknot [14, 24-1. If so, output a pre- constructed efficient suitable 
triangulation of the unknot complement and terminate immediately. 

2. Run the procedure of Hass, Lagarias and Pippenger 1151 Lemma 7.2] to obtain a triangu- 
lation T of K , along with a meridian expressed as a path that follows boundary edges of 

r. 

3. Run the procedure of Jaco and Rubinstein 1211 Proposition 5.15 and Theorem 5.20] to 
convert this into a one-vertex triangulation with no embedded normal 2-spheres. Keep 
track of the location of the meridian as the Jaco- Rubinstein procedure runs. 

4- Apply layerings of extra tetrahedra to alter the boundary edges until the meridian consists 
of a single boundary edge, and output the resulting triangulation. 

Each layering attaches a new tetrahedron along two boundary faces as illustrated in Fig- 
ure^ The number and locations of these layerings are determined by a continued fraction 
calculation, as described by Jaco and Rubinstein 122\ \23\[ . 

new tetrahedron 




boundary 



Figure 4: Layering a new tetrahedron onto the boundary 

The unknot test in step [1] is necessary because the Jaco-Rubinstein procedure in step [3] 
requires K to be both irreducible and boundary irreducible, which is only true for non-trivial 
knots [15] . If K is the unknot, the one-tetrahedron solid torus [21] can be used as a pre- 
contracted solution. 



Theorem 3.2. Algorithm \3. 1\ is correct, i.e., it produces an efficient suitable triangulation of 
K. 

Proof. For the correctness of each sub-procedure we refer the reader to the source papers [14"] [15 ] 
I2T1 [221 [231 124"] . Here we simply prove that the requirements for an efficient suitable triangulation 
are satisfied. 

The Jaco-Rubinstein procedure in step [3] gives a triangulation with no embedded normal 
2-spheres and just one vertex. The only missing requirement is that some boundary edge is a 
meridian. The layering process in step [4] fixes this, and does not break the other requirements: 

• Each layering preserves the number of vertices of the triangulation. 

• The layering process does not introduce any new embedded normal 2-spheres. Each time 
we layer a new tetrahedron onto the boundary, every normal triangle or quadrilateral 
in this new tetrahedron meets the boundary of the new triangulation. Therefore any 
embedded normal 2-sphere in the new triangulation cannot use these new normal discs, 
and must have been an embedded normal 2-sphere in the old triangulation also. □ 



Steps Q] and [3] of the previous algorithm are slow; although step Q] can be avoided (e.g., by 
using prior information that the input knot is non-trivial), the Jaco- Rubinstein procedure of 
step [3] remains too inefficient to use with all but the simplest knot complements. 

We therefore offer an alternative, heuristic algorithm that uses only fast (small polynomial 
time) operations. The drawback is that this heuristic algorithm might not produce any solution 
at all; however, experience shows this to be a rare occurrence, as discussed below. 

Algorithm 3.3. Given a knot diagram, describing the knot K C S 3 , the following heuristic 
procedure will either output a suitable triangulation of K or terminate with no output at all. 

1. Run the procedure of Hass, Lagarias and Pippenger \15\, Lemma 1.2} to obtain a triangu- 
lation T of K , along with a meridian expressed as a path that follows boundary edges of 
T. 

2. Simplify the triangulation using fast local operations (such as edge collapses, Pachner 
moves, book closing moves and related operations 14- 91) to reduce the number of tetrahedra 
and the number of boundary faces as far as possible. 

3. Test whether the resulting triangulation T' is suitable. If so then output T' , and otherwise 
terminate with no output. 

The choice of local simplification operations in step[2]is not important; see [HE] for details. 
Although the suitability test in step [3] requires 3-sphere recognition (which runs in worst-case 
exponential time), we can sidestep this by tracking the location of the meridian throughout 
step El avoiding the need to verify the meridian condition (and test for 3-spheres) in step [3j 

In Section [7] we observe that this heuristic algorithm outputs a suitable triangulation for 
all 2977 non-trivial prime knots with < 12 crossings, showing it to be extremely effective in 
practice. 

We finish this section with some useful properties of suitable triangulations. 

Lemma 3.4. Let T be an efficient suitable triangulation of a knot complement. Then every 
closed normal surface embedded in T has Euler characteristic < 0. 

Proof. The efficiency criterion ensures that there are no embedded normal 2-spheres. The only 
other closed surface of positive Euler characteristic is the projective plane, which does not embed 
in K 3 and so cannot embed in any knot complement. □ 

Lemma 3.5. Let K be any knot, let T be any suitable triangulation of its complement with 
meridional boundary edge m, and let S be any normal surface in T ■ Then S is a spanning 
surface for K if and only if S has no closed components and S meets edge m in precisely one 
point. 

Proof. Let dT denote the two-triangle torus on the boundary of T, and let dS denote the 
collection of boundary curves of S. Each curve of dS is a normal curve on dT] that is, a union 
of arcs between distinct edges of dT as illustrated in Figure [SJ 

If S is spanning then S is connected and has no closed components; moreover, dS consists 
of a single normal curve c whose algebraic intersection number with m is ±1. It is a property 
of normal curves on a two-triangle torus that every such curve meets m in precisely one point 
[23]. 

Conversely, suppose that S has no closed components and meets m in precisely one point. 
Then S has some boundary curve c whose algebraic intersection number with m is ±1. Because 






Figure 5: Examples of normal curves on a two-triangle torus 

the boundary curves of S are disjoint and dT is a torus, any other boundary curve d of S must 
be parallel to c or trivial in dT, both of which would generate additional intersections with m 
(note that the only trivial normal curve on a two-triangle torus cuts each edge twice). Therefore 
c is the only boundary curve of S. 

Since S has no closed components it follows that S is a connected surface with boundary c, 
and since c has algebraic intersection number ±1 with m it follows that S is spanning. □ 

4 A Hilbert basis algorithm 

Our first algorithm for computing the crosscap number C{K) follows a common pattern for 
topological algorithms: it is based on the enumeration of fundamental normal surfaces. 

It is worth revisiting the algorithm for the orientable counterpart of C(K), the knot genus 
g(K), as described by Hass, Lagarias and Pippenger [15]. Their algorithm is based on the 
following result: 

Theorem 4.1 (Hass, Lagarias and Pippenger, 1999). Let K be any knot, and T be any trian- 
gulation of its complement. Then there is a fundamental normal orientable spanning surface of 
genus g{K). 

The algorithm for computing g{K) is then to enumerate all fundamental normal surfaces in 
T, and to observe the smallest genus orientable spanning surface that appears. 

For computing crosscap number, things are less straightforward: it is not even known whether 
there must be a normal non-orientable spanning surface of non-orientable genus C{K), let alone 
a fundamental normal surface. The arguments of Hass, Lagarias and Pippenger use the fact 
that any minimal genus orientable spanning surface is essentiaQ however, these arguments do 
not translate to the non-orientable case since there are knots for which every minimal genus 
non-orientable spanning surface is non-essential [H 119] . 

Our solution is twofold: we work with efficient suitable triangulations, which allow us to 
obtain precise results in many cases, and for those cases that remain we use Clark's inequality 
to reduce the solution to one of two possible values. 

Lemma 4.2. Let K be any non-trivial knot, and T be any suitable triangulation of its comple- 
ment. Then either there is a normal non-orientable spanning surface of non-orientable genus 
C(K), or else C{K) = 2g(K) + 1. 

Proof. Let S be any non-orientable spanning surface of non-orientable genus C(K), and let m be 
the meridional boundary edge of T. Since the boundary of S has algebraic intersection number 
±1 with ra, we can isotope S so that the boundary of S meets m in precisely one point (a simple 
operation on the two-triangle torus boundary of T) ■ 



4 That is, both incompressible and boundary incompressible. Further details are not required here. 



We now follow the standard normalisation procedure that converts an arbitrary properly 
embedded surface into an embedded normal surface (possibly with different topology). We do 
not reiterate the details of this procedure here; for full details the reader is referred to a standard 
reference such as [20] or [25]. Instead we highlight some of its key aspects: 

• Most steps in the procedure are isotopies, which preserve the fact that S is a non-orientable 
spanning surface of genus C(K). However, three types of step can alter the topology of S: 

— internal compressions, which involve surgery on a disc that is bounded by a curve on 
S, as illustrated in Figure 6(a)| 



boundary compressions, which involve surgery on a disc that is bounded by an arc on 
S and an arc on the boundary of T, as illustrated in Figure |6(b)[ 

deletion of trivial components, where we remove components of S that are trivial 
spheres (bounding a ball in K) or trivial discs (parallel into the boundary of K) . 





(a) An internal compression 
^ — boundary of T 

disc S 

^±zz^ — ► J^ ^ 

(b) A boundary compression 

Figure 6: Examples of normalisation moves that can alter the topology of S 

• For any edge e of T, no step in the procedure will ever increase the number of intersec- 
tions between the surface and e. Moreover, each internal compression and each boundary 
compression can be performed in a manner that never reduces this number of intersections 
either. 

We now analyse how internal compressions and boundary compressions can affect our non- 
orientable spanning surface S. At each stage, we assume that S is a non-orientable spanning 
surface of minimum genus C(K). 

• For internal compressions, there are two cases to consider: 

— If the boundary of the compression disc is a separating curve in S, then the com- 
pression splits S into two disjoint pieces Si and Sj, where S is the connected sum 
Si^S-j. Without loss of generality, we assume that the boundary curve of S stays 
with Si, and so S\ is also a spanning surface and S% is a closed surface. 
Since S = Si#S2, at least one of Si and S2 must be non-orientable, and since no 
closed non-orientable surface can embed in a knot complement, this non-orientable 
piece must be Si- Moreover, a simple Euler characteristic calculation gives x(S) = 
x(Si) + XOS2) — 2, and since x(&2) — 2 it follows that Si (like S) must have the 
smallest possible non-orientable genus C(K). Therefore we can simply delete S2, 



replace 5 with the non-orientable spanning surface 5*1 of genus C(K), and continue 
the process. 

— If the boundary of the compression disc is not separating in 5, then the resulting 
surface is some spanning surface U for which x(U) — x(S) + 2. U cannot be non- 
orientable, since its non-orientable genus would be less than the minimum C(K). 
Therefore U is an orientable spanning surface with orientable genus |(1 — x{U)) = 
|(1 - x(S) - 2) = \{C(K) - 2). This gives g(K) < \{C{K) - 2) in contradiction to 
Clark's inequality (Theorem 12. ip , and so this case can never occur. 

• For boundary compressions we first note that, since boundary compressions do not change 
the number of intersections with the meridional edge m, the resulting (possibly discon- 
nected) surface still cuts m precisely once. 

— If the compression separates 5 into two disjoint surfaces S± and 52, we may assume 
that Si cuts edge m precisely once and 52 does not cut m at all. Therefore the 
boundary of 52 must be either a meridian or a trivial curve on the boundary torus; 
either way, it can be filled with a disc in S 3 \K to yield a closed surface in 5 3 . Since no 
closed non-orientable surface can embed in 5 3 , it follows that 52 must be orientable. 
Therefore the piece 5i is non-orientable. This time the Euler characteristic argument 
gives x(5) = x(Si) + xifo) — 1 with x(5 , a) < 1- Therefore the piece Si is again a 
non-orientable spanning surface with genus C(K), and we simply delete 52, replace 
5 with Si and continue. 

— If the compression does not separate 5, the result is a spanning surface U with x(^0 = 
x(S) + 1. In this case the normalisation process might fail. As before, U cannot have 
non-orientable genus less than C(K); therefore U is an orientable spanning surface 
with orientable genus |(1 - x(U)) = \{C{K) - 1). This gives g{K) < \{C{K) - 1), 
and combined with Clark's inequality we obtain C(K) — 2g(K) + 1 precisely. 

If we follow the arguments above, deletion of trivial components is never required: we delete 
unwanted extra components as they appear during internal and boundary compressions, and 
the spanning surface itself is never a trivial sphere nor a trivial disc. 

In conclusion, either the normalisation procedure yields an embedded normal non-orientable 
spanning surface of genus C(K), or else we have a situation in which C(K) = 2g(K) + 1. □ 

We can now move from normal surfaces to fundamental normal surfaces, giving us the main 
theorem of this section. 

Theorem 4.3. Let K be any non-trivial knot, and T be any efficient suitable triangulation of 
its complement. Then either there is a fundamental normal non-orientable spanning surface of 
non-orientable genus C(K), or else C(K) s {2g(K), 2g(K) + 1}. 

Proof. If there is no normal non-orientable spanning surface of genus C(K) then C(K) = 
2g(K) + 1, by Lemma 14.21 Assume then that there is some normal non-orientable spanning 
surface of genus C(K), and let 5 be such a surface containing the fewest possible normal discs. 
Let m denote the meridional boundary edge of T ■ 

If 5 is not fundamental then 5 = U + V for some non-empty normal surfaces U and V. 
By Lemma 13.51 5 cuts the edge m in precisely one point. Since edge intersections are additive 
under normal sun£J, we may assume without loss of generality that U cuts m precisely once 



5 This refers to absolute numbers of intersections, not algebraic intersection number. 
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and V does not cut in at all. We may also assume that U is connected (since any components 
that do not meet m can be moved over to V ) , whereby Lemma 13.51 shows that U is a spanning 
surface also. 

The surface V is a union of closed components and/or bounded components. Every boundary 
curve of V is disjoint from m and is therefore a meridian (since any other normal boundary 
curve, even the trivial curve, must cut m at least once). This means that no bounded component 
of V can be a disc (because the meridian has non-trivial homology in K, and so no disc in K can 
have a meridian as boundary) . Therefore every bounded component of V has non-positive Euler 
characteristic; combined with Lemma [3.41 this gives us x(Y) <■ 0- Since Euler characteristic is 
additive under normal sum we have x{S) = x(U) + x(V), and therefore x{U) > x(S)- 

If U is non-orientable, it follows that U is a normal non-orientable spanning surface of genus 
at most C(K). By minimality of C(K) this genus is precisely C(K); moreover, since S = U + V 
we see that U contains fewer normal discs than S. This contradicts our initial choice of 5*. 

Therefore U is orientable, and has genus < i(l — x{S)) = \C(K). Since U is spanning this 
gives 2g{K) < C(K), and by Clark's inequality it follows that C{K) G {2g(K), 2g(K) + 1}. □ 

Combining Theorems 14.11 and 14.31 gives our first algorithm for computing crosscap number. 

Algorithm 4.4. Given a knot diagram describing the knot K C S 3 , the following procedure 
will either (i) output the crosscap number C(K), or (ii) output a pair of integers one of which 
isC(K). 



1. Construct an efficient suitable triangulation T of the complement K using Alaorithm \3.1\ 

2. Enumerate the set T of all fundamental normal surfaces in T '. 

3. For each surface S G T , test whether S is a spanning surface for K and whether S is 
orientable or non-orientable. Let g n be the minimum non-orientable genus amongst all 
non-orientable spanning surfaces in T , or oo if no non-orientable spanning surfaces are 
found. Let g be the minimum orientable genus amongst all orientable spanning surfaces 
in J- '. 

4- If go = then output 0. If g n < 2g then output g n . Otherwise output {2g 0l 2g + 1}. 

As outlined in [T3] , the enumeration in step [5] involves a high-dimensional Hilbert basis 
computation, and can be performed using standard software such as Normaliz [3] . In step [3l 
we can test for orientability and genus by reconstructing the normal surface, and we can test 
whether S is a spanning surface for K using Lemma [ 



Theorem 4.5. Algorithm ^ -4\ is correct, i.e., the crosscap number C(K) is one of the solutions 
that it outputs in stepU\ 

Proof. By Theorem 14. 11 g = g(K). If K is the unknot then g = and the algorithm outputs 
the correct value C(K) = 0. Assume then that K is non-trivial, which means that g Q > 0. 

If some non-orientable spanning surface of genus C(K) does appear in the set T then we also 
have g n = C(K). If g n < 2g then we output the correct value g n — C(K). Otherwise Clark's 
inequality (Theorem I2.1J) gives C(K) = 2g + 1, and so the correct value of C(K) appears in 
the output set {2g , 2g + 1}. 

If no non-orientable spanning surface of genus C(K) appears in T then g n > C(K), and by 
Theorem l4~3l we also have C{K) G {2g(K), 2g{K) + 1} = {2g , 2g + 1}. In this case 2g < g n , 
and again our output set {2g , 2g + 1} contains the correct value of C(K). □ 
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5 An exact integer programming algorithm 

The bottleneck in Algorithm 14.41 is the Hilbert basis computation (the enumeration of all fun- 
damental normal surfaces), which remains intractable for all but the simplest knots. For this 
reason we seek alternate algorithms that do not require an enumeration of surfaces. Instead 
we attempt to locate a minimal genus non-orientable spanning surface directly using integer 
programming. 

We begin this section with several results that allow us to formulate topological constraints 
numerically. We finish with Algorithm 15 .41 which gives the full procedure for computing C{K). 

The first lemma is due to Hass, Lagarias and Pippenger [T5], and allows us to place upper 
bounds on the coordinates of fundamental normal surfaces. 

Lemma 5.1 (Hass, Lagarias and Pippenger, 1999). Let T be a 3-manifold triangulation with 
n tetrahedra, and let S be a fundamental normal surface in T . Then every coordinate of the 
vector representation v(5 l ) is at most n ■ 2 +2 . 

This bound plays an important role in our algorithm, and it is desirable to improve it. Tech- 
niques based on linear programming offer potential; for instance, in [8] the authors significantly 
reduce this bound for related problemsO We do not pursue these matters further here. 

Our next result is the well-known observation that Euler characteristic is linear on the 
normal surface solution cone. Our proof includes details on how such a linear function can be 
constructed. 

Lemma 5.2. Let T be an n-tetrahedron 3-manifold triangulation. Then there is a linear func- 
tion x : K 7n — y R such that, for any normal surface S in T, x(v(S*)) is the Euler characteristic 
ofS. 

Proof. Let x G R 7n . We define x(x) = V - E - F, where: 

• F is the sum of all 4n triangular coordinates and all in quadrilateral coordinates of x. 

• E is defined as a sum over edges of normal discs: if e is one of the three edges of a normal 
triangle or one of the four edges of a normal quadrilateral, then e contributes +1 to this 
sum if e lies on the boundary of T, or +1/2 if e is internal to T ■ 

• V is defined as a sum over vertices of normal discs: if v is one of the three vertices of a 
normal triangle or one of the four vertices of a normal quadrilateral, then v contributes 
+l/d to this sum, where d is the degree of the edge of T that v lies within. 

It is clear that V, E and F can all be expressed as sums of normal coordinates, and so x is 
linear on M. 7n . Moreover: each vertex of degree d in a normal surface S contributes a total of 
+ l/d x d = +1 to V (since it meets d normal discs, and it must lie within an edge of 7~ of 
degree d); each boundary edge of S contributes +1 to E; each internal edge of S contributes 
+1/2 x 2 = +1 to E (since it meets two normal discs); and each triangular or quadrilateral face 
of S contributes +1 to F. Therefore x(v(S')) = V — E + F is indeed the Euler characteristic of 
S. □ 

It should be noted that there are many linear functions \ '■ ^ 7n —> K with this property; the 
formulation above was chosen for its simple proof. There are sparser formulations that may be 
preferable for computation; again we do not pursue this matter here. 



6 The paper [8] considers vertex (not fundamental) normal surfaces, and works in a different coordinate system. 
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Our final lemma allows us to arithmetically express the condition that a given normal surface 
is a spanning surface for our knot. 

Lemma 5.3. Let K be any knot, and let T be any suitable triangulation of its complement with 
meridional boundary edge m. Let A$ be any tetrahedron of 1~ that contains the boundary edge 
ra, and let ti ia > ti,&, qi.c and qi t d be the coordinates describing the four normal disc types in Aj 
that touch m, as illustrated in FigureVK 

For any normal surface S in T, S is a spanning surface for K if and only if S has no closed 
components and the vector representation of S satisfies tt, a + ti,b + Qi,c + Qi,d = 1- We call this 
equation the spanning equation for T ■ 







Figure 7: The four normal disc types in Aj that touch the meridional boundary edge m 

Proof. The sum ti, a + ij,& + qi, c + <k,d counts the number of times that the surface S cuts the 
meridional boundary edge m (this is true regardless of which tetrahedron A, we chose) . The 
result then follows immediately from Lemma 13.51 □ 



Algorithm 5.4. Given a knot diagram describing the knot K C S 3 , the following procedure 
will either (i) output the crosscap number C{K), or (ii) output a pair of integers one of which 
isC{K). 



1. Construct an efficient suitable triangulation T of the complement K using Alaorithm \3.1\ 
Let n be the number of tetrahedra in T ■ 

2. Optimise the following integer program using an exact arithmetic integer programming 
solver: 

(i) define the In non-negative integer variables tjj for 1 < % < n, 1 < j < 4 and qij for 
1 < i < n, 1 < j < 3; 

(ii) define the 3n binary variables bij for 1 < i < n, 1 < j ' < 3; 

(Hi) add constraints for the matching equations and the spanning equation (Lemma \5.3\) : 

(iv) add constraints n ■ 2 7n+2 • bij > q%j for 1 < i < n, 1 < j < 3; 

(v) add constraints bi\ + bifl + &j,3 < 1 for 1 < i < n; 

(vi) maximise the Euler characteristic function (Lemma \5.2\) . 



3. Let x be an optimal solution, and let S be the corresponding normal surface with any closed 
components removed. If S is non-orientable then output 1 — x(S'). Otherwise output the 
patr{l-x(S),2- X (S)}. 



13 



It is crucial that the integer programming solver be based on exact integer arithmetic; see 
[TJ [T2] for examples of such tools. With traditional floating-point solvers, round-off errors can 
creep in, especially when working with coefficients as large as n ■ 2 7 "+ 2 . Such round-off errors 
can result in a sub-optimal solution, or even a solution x that does not represent a spanning 
surface at all. 

Theorem 5.5. Algorithm \5.4\ is well-defined (i.e., the integer program in step\^ is not un- 
bounded), and is correct (i.e., the crosscap number C(K) is one of the solutions output in 
stepWj). 

Proof. The integer variables ti t j and qt.j correspond to the usual triangle and quadrilateral 
coordinates. Each binary variable bij has the following effect under condition (iv): if bij = 
then it forces qij = 0, and if bij = 1 then we can have any q. L ,j in the range < qij < n ■ 2 7 "+ 2 . 

First consider any feasible solution to the integer program (optimal or otherwise). Because 
any such solution satisfies all constraints from step [5J we know that the variables tij and q^j are 
non- negative integers, satisfy the matching equations and the spanning equation (condition (iii)), 
and satisfy the quadrilateral constraints (condition (v)). By Theorem 12.21 and Lemma [5.31 it 
follows that any feasible solution to our integer program describes a normal surface in T, which 
is the disjoint union of a spanning surface and zero or more additional closed components. 

If the integer program is unbounded then such a union can have arbitrarily large Euler 
characteristic. Because any spanning surface has Euler characteristic < 1, it follows that T 
contains some closed normal surface of positive Euler characteristic. However, this is impossible 
by Lemma 13.41 and so the integer program is bounded and Algorithm 15 .41 is well-defined. 

We now turn to the proof of correctness. By the argument above, the surface S obtained 
in step [3] (with closed components removed) must be a spanning surface. If S is non-orientable 
then its non-orientable genus is 1 — x(S), and we have C(K) < 1 — x{S)- If S is orientablc 
then its orientable genus is |(1 — x(S)) > g{K), and by Clark's inequality we have C(K) < 
2g(K) + l<2- X (S). 

Having bounded C(K) from above, we now bound it from below. Let S be the smallest- 
genus orientable normal surface whose normal coordinates are all at most n ■ 2 7n+2 , and let 
S n be the smallest-genus non-orientable normal surface whose normal coordinates are at most 
n ■ 2 7n+2 . By Theorem 14.11 and Lemma l5.1[ the surface S exists and has genus g(K). By 
Theorem 14.31 and Lemma 15. 1[ we have one of two cases for S n '■ either (a) S n exists and has 
non-orientable genus C(K), or (b) C(K) > 2g(K) (in which case S n might or might not exist). 

Since S and S n are normal spanning surfaces with coordinates < n ■ 2 7n+2 , they satisfy all 
of the constraints in step [5] of the algorithm (here each bij = or 1 according to whether the 
corresponding q, t j is zero or non-zero). If S" is the normal surface corresponding to the optimal 
solution x, it follows that x(So) < x(S'), and that x(^n) < x(S') if S n exists. As before, T 
contains no closed normal surfaces of positive Euler characteristic, and so after removing any 
closed components of S' we obtain x(So) < x(S), and x(S'„) < x{S) if S n exists. 

We can now piece our various results together. In case (a) above where S n exists with 
non-orientable genus C(K), we have C(K) = 1 — x(5„) > 1 — x(S'). In case (b) we have 
C(K) > 2g{K) = l-x(So) > l-x(S). Either way we obtain the lower bound 1 ~x(S) < C{K). 

In conclusion: if S is orientable then 1 — x(S) < C(K) < 2 — x(S), and if S is non-orientable 
then 1 — x(S) < C(K) < 1 — x(S'). Therefore the output of our algorithm is correct. □ 
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6 A limited precision integer programming algorithm 

The exact integer programming algorithm in the previous section avoids an expensive Hilbert 
basis enumeration, but it has its own drawbacks. Exact integer programming solvers are rarer 
and less well-developed than their limited precision floating-point cousins, and the exponentially 
large constraint coefficients n ■ 2 7n+2 can have a crippling effect on performance. Moreover, 
constructing an efficient suitable triangulation remains slow, as discussed in Section G3 

Our final algorithm avoids these performance problems: we allow floating-point solvers with 
limited precision, and we replace each large coefficient n ■ 2 +2 in our integer program with the 
arbitrarily chosen small coefficient fOOOO. We also drop the efficiency requirement and allow 
just suitable triangulations, which are faster to construct. 

These concessions bring about a loss of information and precision from our solution. In 
response, we add tests to ensure that the solution to our integer program is valid (i.e., represents a 
spanning surface), and we interpret the final output value as just an upper bound on C(K) (since 
we cannot be sure that a better solution was inadvertently missed). The complete algorithm is 
as follows. 

Algorithm 6.1. Given a knot diagram describing the knot K C S 3 , the following procedure will 
output an upper bound U for which C(K) < U . 

1. Construct a suitable triangulation T of the complement K using the fast but heuristic-based 
Algorithm \3.3i If this algorithm produces no triangulation then output oo and terminate 
immediately. Otherwise let n be the number of tetrahedra in T . 

2. Optimise the following integer program using a fast solver based on floating-point arith- 
metic: 

(i) define the In non-negative integer variables t% t j for 1 < i < n, 1 < J < 4 and qij for 
1 < i < n, 1 < j < 3; 

(ii) define the 3n binary variables bi.j for 1 < i < n, I < j < 3; 

(Hi) add constraints for the matching equations and the spanning equation (Lemma \5.3\) ; 
(iv) add constraints 10000 • bij > qij for 1 < i < n, 1 < j < 3; 

(v) add constraints 6^i + 6^,2 + 6i,3 < 1 for 1 < i < n; 
(vi) maximise the Euler characteristic function (Lemma \5.2\) . 

3. If this integer program is unbounded then output oo and terminate. 

4-. Let x be an optimal solution. Using exact integer arithmetic, test whether x satisfies the 
constraints laid out in step // not then output ex) and terminate. 

5. Let S be the normal surface described by x with any closed components removed. Output 
1 — x{S) if S is non-orientable, or 2 — x(S) if & * s orientable. 

Theorem 6.2. Alaorithm \6.1\ is correct, i.e., C(K) is less than or equal to the output value. 



Proof. This is a much simpler variant of the proof for Algorithm 15.41 If we output oo then the 
algorithm is clearly correct. Otherwise we have a solution x that satisfies the constraints of 
stepd] (as shown by the test in stepQJ, but with no guarantee that this solution is optimal. 
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As in Algorithm 15.41 because x satisfies the constraints from step [H the surface S (with 
closed components removed) must be a spanning surface for K. If S is non-orientable then its 
non-orientable genus is 1 — x(S'), and so C{K) < 1 — x(S)- If S is orientable then its orientablc 
genus is ~(1-%(S)) > g(K), and by Clark's inequality we have C(K) < 2g(k)+l < 2-%(5). D 

7 Discussion and computational results 

As discussed in the introduction, Algorithms 14.41 and 15.41 remain too slow for practical use for 
all but the simplest knots. Algorithm 14. 41 requires the enumeration of a Hilbert basis for a high- 
dimensional polyhedral cone, an extremely expensive procedure. Algorithm 15.41 requires exact 
arithmetic integer programming in problems with extremely large coefficients (such as n-2 7n+2 ), 
which remains intractable with currently available tools. 

Despite this, both algorithms have certain benefits. Unlike our final algorithm, they produce 
no more than two possible solutions; moreover, further analysis shows that if C(K) < 2g(K) 
then both Algorithms 14.41 and 15.41 guarantee a unique solution. They also rely on substantially 
different underlying computational problems (Hilbert basis enumeration versus exact integer 
programming), both of which have supporting software that continues to enjoy significant ad- 
vances in efficiency [TJ [31 [T2] ; in this sense we are able to "hedge our bets" . A further avenue 
for improving the efficiency of Algorithm 15.41 is to lower the coordinate bounds in Lemma I5TT1 

An important observation is that none of the algorithms in this paper are able to give a 
unique solution in cases where Clark's inequality is tight, i.e., C(K) = 2g(K) + 1. 

In contrast to the first two algorithms, the limited-precision Algorithm 16. H is fast and effec- 
tive. By combining the upper bounds from this algorithm with lower bounds from existing knot 
tables, we are able to make significant improvements to these tables. 

Specifically: the Knotlnfo project contains a rich body of invariant data for all 2977 prime 
non-trivial knots with < 12 crossings. Of these knots, 289 have known crosscap numbers, and 
the remaining 2688 knots are listed with best-known upper and lower bounds. 

For each of these knots, we begin with the corresponding triangulation of the complement 
from the SnapPy census [T3] , use Regina [U [7] to convert it into a triangulation with boundary 
faces, and then run Algorithm 16.11 to obtain a bound on the crosscap number. For the crucial 
optimisation step we use IBM's ILOG CPLEX package (version 12.2). The triangulations range 
from 5 to 50 tetrahedra in size (with an average of 33.27), and the total running time over all 
triangulations is roughly 4.5 hours on a quad-core 2.93 GHz Intel Core i7 CPU. 

Our first observation is that Algorithm 16. II never outputs oo. That is: for every knot in the 
tables, the heuristic Algorithm 13.31 produces a suitable triangulation, the integer program is not 
unbounded (giving strong evidence that these triangulations may in fact be efficient suitable 
triangulations), and the optimal solution to the integer program satisfies all of the necessary 
constraints to produce a valid spanning surface. 

For 27 of the 2688 knots with unknown crossing number, the upper bound produced by 
Algorithm 16.11 is equal to the lower bound listed in the Knotlnfo tables; as a result we can 
identify the crosscap number precisely. These 27 knots and their crosscap numbers are listed in 
Table [T] For another 747 knots, our upper bound improves upon the Knotlnfo upper bound: 
on average we reduce the number of possible solutions from 4.18 to 2.86. Detailed results from 
all of these computations, including a Regina data file listing the final spanning surfaces, can be 



found at http://www.maths.uq.edu.au/~bab/code/. 



To conclude, we note that the integer programming framework given in Algorithms 15.41 and 
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Knotlnfo 


Dowker-Thistlethwaite 


Genus 


Previous bounds 


New value 


name 


name 




on C(K) 


of C(K) 


820 


8m 


2 


[2,4] 


2 


10ia8 


10ni 5 


3 


[2 


4] 


2 


10x26 


IO1117 


3 


[2 


4] 


2 


IO139 


10n27 


4 


[2 


3] 


2 


10X40 


10n 2 g 


2 


[2 


4] 


2 


10x42 


10n 30 


3 


[2 


4] 


2 


10X45 


lOmx 


2 


[2 


4] 


2 


10l61 


10n 3 i 


3 


[2 


5] 


2 


Hnio2 


Hnio2 


2 


[2 


4] 


2 


1111X04 


II11104 


4 


[2 


4] 


2 


H11135 


Ilni3s 


3 


[2 


5] 


2 


12n 121 


12noi2x 


2 


[2 


4] 


2 


12ll0233 


12no233 


4 


[2 


4] 


2 


12ll0235 


12no235 


4 


[2 


4] 


2 


12ll0242 


12no242 


5 


[2 


3] 


2 


12no404 


12no404 


2 


[2 


4] 


2 


12ll0474 


12no474 


4 


[2 


4] 


2 


12ll0475 


12no475 


3 


[2 


4] 


2 


12ll0522 


12no522 


3 


[2 


4] 


2 


12n 575 


12n0575 


4 


[2 


4] 


2 


12no58i 


12n 58i 


3 


[2 


4] 


2 


12n 582 


12no582 


2 


[2 


4] 


2 


12no59i 


12no59i 


4 


[2 


5] 


2 


12no72i 


12no72i 


4 


[2 


4] 


2 


12ll0725 


12no725 


5 


[2 


3] 


2 


12ll0749 


12no749 


3 


[2 


5] 


2 


12no85i 


12n 85i 


3 


[2 


5] 


2 



Table 1: The 27 knots with newly-computed crosscap numbers 



6.11 extends beyond the specific problem of computing crosscap numbers — it can be used as a 
general framework for locating normal surfaces with various properties. For instance, by optimis- 
ing Euler characteristic this integer programming framework can be used to test for 0-efhciency 
[2"T] . compute connected sum decompositions [21] . and perform 3-sphere recognition [5]. Given 
current advances in exact integer programming solvers, this potential for applying optimisation 
techniques to problems in low-dimensional topology is only beginning to be explored. 
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